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The effect of layer-by-layer heterostructuring and epitaxial strain on lattice instabilities and re- 
lated ferroelectric properties is investigated from first principles for the [001]-(PbTiO3)i(PbZrO3)i 
superlattice and pure PbTiC>3 on a cubic substrate. The results for the superlattice show an enhance- 
ment of the stability of the monoclinic r-phase with respect to pure PbTiC>3. Analysis of the lattice 
instabilities of the relaxed centrosymmetric reference structure computed within density functional 
perturbation theory suggests that this results from the presence of two unstable zone-center modes, 
one confined in the PbTiC>3 layer and one in the PbZrC>3 layer, which produce in-plane and normal 
components of the polarization, respectively. The zero-temperature dielectric response is computed 
and shown to be enhanced not only near the phase boundaries, but throughout the r-phase. Analysis 
of the analogous calculation for pure PbTiC>3 is consistent with this interpretation, and suggests 
useful approaches to engineering the dielectric properties of artificially structured perovskite oxides. 
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I. INTRODUCTION 

Perovskite oxide superlattices are expected to have 
properties distinct from those of the bulk constituents. 
Due to lattice mismatch, constituent layers below the 
critical thickness generally are in highly strained states. 
As can be confirmed by first-principles studies of bulk 
materials, the lattice instabilities and properties of bulk 
perovskites are extremely sensitive to the strain state. 
The degree of strain can be further tuned by varying the 
lattice constant of the substrate on which the superlattice 
is grown. 

This discussion is of practical interest as it has been 
shown that perovskite oxide superlattices can indeed be 
grown with atomic scale precisionii2*£, opening up a 
highly varied family of artificially structured materials 
for investigation. Thus, first-principles prediction of su- 
perlattice properties can directly lead to productive inter- 
actions with experimentalists and progress in materials 
design^. 

In this work, we consider the simplest possible (and 
as yet hypothetical) [001]-(PbTiO 3 ) m (PbZrO 3 )„ super- 
lattice with m=l and n=l. As shown in Fig.^ this con- 
sists of single unit cell layers of PbZr03 (PZ) and PbTiC>3 
(PT) alternating along [001]. The equilibrium cubic lat- 
tice constant of bulk PZ is significantly greater than that 
of PT, with the computed values being ao.pz=7.77 a.u. 
and ao, pt=7. 37 a.u., leading to a 5% mismatch (the cor- 
responding experimental values are 7.81 a.u., 7.50 a.u, 
and 4%). At their equilibrium lattice constants, both cu- 
bic PT and cubic PZ have an unstable three-fold degen- 
erate polar phonon mode at the zone center. The lattice 
mismatch means that in the superlattice, PZ is under 
compressive in-plane stress, while PT is under tensile in- 
plane stress. If the structure is allowed to relax without 
breaking the central mirror plane symmetry (i.e. in space 
group P4/mmm), the in-plane lattice constant takes an 
intermediate value of 7.57 a.u.. The PZ unit cell elon- 




FIG. 1: The ideal double-perovskite 10-atom unit cell of the 
[001]-(PbTiO 3 )i(PbZrO 3 )i superlattice. The lead, oxygen, 
titanium and zirconium atoms are indicated by small solid 
circles, light grey circles, dark grey circles and open circles, 
respectively. The oxygen octahedra and the ideal perovskite 
5-atom cells are outlined. The pseudoperovskite c/a is defined 
as half the unit cell c/a. Note that this is the same structure 
considered as a ordered realization of PZT in Ref. 0. 



gates, giving a local tetragonal unit cell with c/a > 1, 
and conversely, the PT unit cell will shorten, resulting 
in a local tetragonal unit cell with c/a < 1. This low- 
ering of unit cell symmetry leads to a splitting of the 
three-fold degenerate unstable mode in each layer, with 
the lower frequency mode(s) along the long direction(s). 
Thus, on qualitative grounds, one expects at least two 
unstable T modes of different symmetry, one associated 
with polarization along the normal in the PZ layer, and 
a second, two-fold degenerate, mode producing polariza- 
tion in the plane of the PT layer. DFPT calculations 
of the phonon dispersion of the relaxed centrosymmetric 
superlattice show two unstable zone-center polar modes, 
one with in-plane polarization confined in the PT layer 
and the other with polarization along the normal confined 
in the PZ layer, quantitatively confirming this picture^. 
The energy of the system can be lowered by either of 
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these modes alone; the actual ground state might involve 
choosing one, or the other, or coupling both together. 

First principles ground-state-structure calculations 
have previously been performed for this system, regard- 
ing the 1:1 [001] ordered PZT 50/50 supercell as an ap- 
proximate realization of the PZT solid solution 7 .. Under 
the constraint of experimental volume, c/a = 1.035 and 
a = (3 = 90 degrees, it was found that the relaxed struc- 
ture has monoclinic Cm symmetry, with internal coor- 
dinates in reasonable agreement with those experimen- 
tally determined for the monoclinic PZT phased. Con- 
necting this result to our discussion in terms of lattice 
instabilities, we see that the ground state is in fact ob- 
tained by freezing in the two modes together, giving a 
total polarization P = (P,P,P Z ) rotated away from the 
c-axis. Piezoelectric coefficients computed for this struc- 
ture in Ref. have very large values. This is attributed 
to the low energy for polarization rotation; this mecha- 
nism was first suggested to explain the high piezoelectric 
response of single crystal PMN-PT2iiS and accounts for 
the enhancement not only of the piezoelectric response 
in the monoclinic phase of the superlattice, but also, as 
we will show below, for the enhancement of the dielectric 
response. 

Under an epitaxial constraint imposed by a cubic [001] 
substrate, the in-plane lattice vectors arc equal in length 
and make an angle of 90 degrees, while the system is oth- 
erwise free to relax. Our qualitative picture of the super- 
lattice instabilities suggests that the two modes should 
respond differently to changes in the substrate lattice 
constant. As the substrate lattice constant increases, the 
mode with polarization along the normal should become 
higher in frequency, while the mode with in-plane polar- 
ization should soften. This change in relative stability 
is expected to affect the direction of the polarization in 
the ground state structure, in a manner consistent with 
previous calculations in which a different structural con- 
straint, on c/a, was imposed. 7 

In this paper, we present a quantitative first-principles 
analysis of the effect of varying in-plane strain on the epi- 
taxially constrained [001]-(PbTiO 3 )i(PbZrO 3 )i superlat- 
tice. In section II, we describe the details of the computa- 
tions. In section III, we present results of the optimized 
structural parameters for in-plane lattice constants a\\ 
ranging from well below do, pt, at which the PT unit cell 
is expected to be approximately cubic, to above ao t pz- 
This wide range of strain allows us more fully to inves- 
tigate the effects of strain on structure and properties. 
As we will see, there is a sequence of phases, from the 
tetragonal c-phase at the smallest values of ay, through 
a monoclinic r-phase with increasing ao and finally to 
an orthorhombic aa-phase. The locations of the phase 
boundaries are determined by computing the Hessian ma- 
trix of each optimized structure, comprised of the zone- 
center force-constant matrix, the elastic constants and 
the quadratic-order coupling between zone-center atomic 
displacements and homogeneous strain, and extrapolat- 
ing to zero eigenvalue. Near the phase boundaries and 



throughout the r-phase the dielectric and piezoelectric 
responses are high, the sensitivity to applied fields and 
stresses resulting from the ease with which the polariza- 
tion can be changed through rotation. We relate these 
results to the lattice instabilities of the centrosymmetric 
reference structure obtained from DFPT computation of 
the q = frequencies. The behavior of the superlattice 
is compared with that of pure PT, which is seen to ex- 
hibit a qualitatively different behavior. In section IV, we 
discuss this distinction using a two-mode vs one-mode 
picture, and suggest some ways that these results can be 
used in materials design. 



II. METHOD 

To predict the ground state structure and zone- 
center phonon frequencies, we use density functional 
theory (DFT) and density functional perturbation the- 
ory (DFPT)Ai*i£ within the local density approximation 
(LDA). 13 Calculations in this work have been done us- 
ing the PWscf packaged We use Vanderbilt ultrasoft 
pseudopotentials 15 to describe the interaction between 
ionic cores and valence electrons, and plane- wave basis 
set with kinetic energy cut-off of 35 Ry. The augmenta- 
tion charges, required by the use of ultrasoft pseudopo- 
tentials, are expanded with an energy cutoff of 350 Ry. 
Brillouin zone (BZ) integrations are performed with a 
6x6x3 Monkhorst-Pack mesh. 

The polarization P is computed using the linearized 
expression P Q = YliB^iap u »/3' wnere u i are the atomic 
displacements of the optimized structure relative to the 
centrosymmetric reference structure (described below), 
and Z* are the atomic Born effective charge tensors. The 
values of polarization reported below are computed with 
the averaged effective charge tensors of pure PT and pure 
PZ given in Table II of Ref. 0, which have been shown to 
be a good approximation to those of the superlattice. 

The epitaxial constraint from a cubic substrate is 
treated implicitly by constraining the in-plane lattice 
constant of the system. We consider a range of in-plane 
strains from -6% to +4%. As the reference in-plane lat- 
tice parameter to define the in-plane strain, we use the 
equilibrium in-plane lattice parameter computed for the 
optimized polar structure, ao = 7.73 a.u., for the super- 
lattice and ao = 7.32 a.u. for PT. Thus, the system at 
zero in-plane strain has the minimum elastic energy. 

For each value of the in-plane strain, we perform the 
structural optimization of the superlattice and of pure 
PT in two stages. First, the ideal double-perovskite 
10-atom unit cell of the superlattice is relaxed with- 
out any symmetry breaking beyond that associated with 
the Zr/Ti ordering (tetragonal space group P4/mmm, 
#123). Similarly, the perovskite 5-atom unit cell of PT is 
relaxed in space group PA/mmm. These will be referred 
to as the centrosymmetric reference structures. The zone 
center phonon modes of these structures are computed 
using DFPT. This shows the presence or absence of un- 
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stable modes that provide a guide to energy-lowering dis- 
tortions. 

Then, the symmetry is broken to allow relaxation into 
a lower symmetry phase, which we label following the no- 
tation of Pertsev et al.ii& For the whole range of in-plane 
strain, we consider relaxations into the c-phase (tetrago- 
nal space group PAmm, #99) with polarization Pz, the 
aa-phase (orthorhombic space group Cmm2, #35), with 
polarization P(x + y), and the r-phase (monoclinic space 
group Cm, #8), with polarization P(x+y)+P z z. We also 
investigated the ac-phase (monoclinic space group Pm, 
#6), with polarization Px + P z z, and the a-phase (or- 
thorhombic space group Pmm2, #25), with polarization 
Px; these two phases are found to be energetically unfa- 
vorable and do not appear in the present phase diagram. 
In principle, an equilibrium c-phase can be obtained by 
optimizing structural parameters from a starting r-phase 
or ac-phase (and an equilibrium aa-phase from a starting 
r-phase) though in practice more accurate results can be 
obtained by studying each of the space groups separately. 
Compatible k-point sampling allows us directly to com- 
pare the energies of these structures. Although the angle 
7 between the c axis and the plane can in principle vary, 
in the total energy calculations it is held fixed at 90 de- 
grees to simplify the calculation. The main effect of this 
approximation is to reduce the range of stability of the 
r-phase as a function of in-plane strain; for the systems 
considered here, we describe below how we compute the 
effect of strain through analysis of the Hessian matrix. 

The zone-center dynamical matrices are computed for 
the optimized structures using DFPT. The zone-center 
phonon modes, obtained by diagonalizing the dynam- 
ical matrices, are used according to the formalism of 
Ref. to compute the phonon contribution to the zero- 
temperature constant-strain static dielectric tensor and 
to study its dependence upon epitaxial strain. We also 
used these results for a systematic study of the stabil- 
ity of the various phases as a function of in-plane strain. 
The full Hessian matrix of the second derivatives of the 
energy with respect to lattice-periodic atomic displace- 
ments and homogeneous strain is constructed with DFPT 
for the atomic displacements and finite difference stress 
and force results for the derivatives involving strain. 
Tracking the lowest eigenvalue of the Hessian matrix, we 
can accurately identify the instability point of a given 
phase as the in-plane strain where the lowest eigenvalue 
becomes negative^. This allows a more accurate deter- 
mination of the position of the phase boundary than total 
energy calculations alone, which are limited by computa- 
tional precision, and also allows us to include the effects 
of the coupling to shear strain 774 and 775 in the r-phase, 
equivalent to allowing 7 to relax from 90 degrees. 



III. RESULTS 

This section is organized as follows. In A, we present 
the calculations of the equilibrium structural parameters 
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FIG. 2: Components of the polarization along the [001] and 
[110] directions, and its magnitude |Ptot], as a function of 
in-plane strain for (a) the [001]-(PbTiO 3 )i(PbZrO 3 )i super- 
lattice, and (b) pure PT. 



and polarization as a function of in-plane strain. Calcu- 
lations of zone-center phonon modes, elastic constants, 
and coupling between the atomic displacements and ho- 
mogeneous strain are performed to carry out the stability 
analysis and precisely locate the phase boundaries. In 
B, we present the results of the zone-center phonons in 
the centrosymmetric reference structure as a function of 
in-plane strain and use them to interpret the phase di- 
agram. In C, we present the in-plane-strain dependence 
of the static dielectric tensor. Results in each subsec- 
tion are presented first for the superlattice and then, for 
comparison, for pure PT. 



A. Structural parameters and phase diagram 

In Fig.|2Ia) we show the structural phase transitions of 
the superlattice with varying in-plane strain. The sym- 
metry is evident from the nonzero components of the 
polarization. There is a continuous (second-order) tran- 
sition from the c-phase, stable at large compressive in- 
plane strain, to the r-phase. In the r-phase, the polariza- 
tion has almost constant magnitude and rotates smoothly 
in the (110) plane from z to x + y, followed by a second 
continuous transition from the r-phase to the aa-phase, 
stable at large tensile in-plane strain. The in-plane-strain 
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FIG. 3: The energy as a function of in-plane strain for the c, aa and r phases, for (a) the [001]-(PbTiO3)i(PbZrOs)i superlattice 
and (b) pure PT. The same scales are used for both systems to facilitate comparison. The insets show the tetragonality c/a as 
a function of in-plane strain. The lowest eigenvalue of the Hessian, for (c) the [001]-(PbTiO3)i(PbZrOa)i superlattice and (d) 
pure PT. The circles are for the c-phase and the triangles for the aa-phase. The solid symbols are the eigenvalues when only 
tetragonal strain is allowed, corresponding to the first- principles results in (a) and (b). The open symbols show the change 
when coupling to shear strain is included, resulting in an enhancement of the stability of the r-phase. 



dependence of the energy in each of these three symme- 
tries is shown in Fig. E^a). The a-phase and ac-phase 
energies computed at ay = 7.57 a.u. were found to be 
significantly higher than that of the r-phase and thus are 
not included in further discussion. Accurate positions of 
the second-order c-r and r-aa phase boundaries can be 
best obtained not from intersection points of the total 
energy curves, but by interpolating the dependence of 
the quadratic-order energy around the equilibrium phase 
at each in-plane strain to obtain the in-plane strain at 
which the instability occurs. 

We compute the Hessian matrix (phonons plus tetrag- 
onal strain, the latter coupling being zero for the c- 
phase and the aa-phase) and plot the lowest eigenvalue 
in Fig. |3Jc) . The zero-crossing of the eigenvalues can be 
seen to correspond to the points where the r-phase be- 
comes more stable in Fig. a) ; in the case of the r-aa 
transition the critical in-plane strain is clearly much eas- 
ier to determine accurately in Fig.HJc). To test the effect 
of fixing 7 to 90 degrees, we also compute the full Hes- 
sian matrix, including coupling to shear strains 774 and 
775, and plot the lowest eigenvalue in Fig. EI C ) • It can be 
seen that the effect of the additional strain relaxation is 
to expand the region of stability of the r-phase, the effect 
being slight for the c-r boundary and more significant for 
the r-aa boundary. 

The corresponding results for pure PbTiC>3 present an 



informative contrast to those for the superlattice. In 
Fig.[3b) we show the structural phase transitions of pure 
PT with varying in-plane strain. The system makes a 
weakly first-order direct transition from the c-phase at 
large compressive in-plane strain to the aa-phase at large 
tensile in-plane strain. The intermediate r-phase present 
in the superlattice has been eliminated. There are small 
discontinuities in the c/a ratio (inset of Fig. EJb)) and 
the magnitude of the polarization at the phase transition. 
The optimized energies for each symmetry as a function 
of in-plane strain are shown in Fig.^b). Optimization in 
the monoclinic space group corresponding to the r-phasc 
produces either the c-phase or the aa-phase depending 
on in-plane strain; these energies lie very slightly higher 
than those obtained by optimizing in the tetragonal space 
groups, reflecting the convergence criterion for the re- 
laxation. We can confirm the weakly discontinuous na- 
ture of this transition by computing the Hessian matrix 
(phonons plus tetragonal strain) and plotting the eigen- 
values in Fig.|3|[d). At the transition in-plane strain, each 
of the two phases is marginally stable. To test the effect 
of fixing 7 to 90 degrees, we also compute the full Hes- 
sian matrix, including coupling to shear strains 774 and 
775, and plot the lowest eigenvalue in Fig. E^d). It can 
be seen that the effect of the additional strain relaxation 
introduces a narrow region of stability of the r-phase, 
between in-plane strains of 0.86% and 1.48%. This win- 
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dow of 0.62% is about three times smaller than the T=0 
width of 2% obtained by Pertsev et al.ji 6 - the latter being 
an extrapolation to low temperature of a Landau theory 
fit near the paraelectric-ferroelectric transition tempera- 
tures above 700 K. 



B. Lattice instabilities of the centrosymmetric 
reference structure 

Next, we investigate the relation of the ground state 
structure to the zone-center phonon modes in the cen- 
trosymmetric reference structure. In the unstrained su- 
perlattice (an =7.57 a.u.), there are two unstable, nearly 
degenerate, polar zone-center modes, an E u mode with 
displacements confined in the PT layer, and an A^ u mode 
confined in the PZ layeri&iS In Fig. 0] we see that these 
behave differently with increasing in-plane strain, the PZ- 
like mode being stabilized by increasing in-plane strain, 
while the PT-like mode becomes increasingly unstable. 
This is consistent with the qualitative picture discussed in 
the introduction, with the mode with polarization along 
the normal being raised in frequency by increasing in- 
plane lattice constant, while the in- plane mode softens. 
The ground state structure reflects these opposing trends 
in the strength of the two unstable modes. At the small- 
est in-plane lattice constant, it is energetically favorable 
to freeze in only the A2 U mode, producing the c-phase. 
As the in-plane lattice constant increases, the in-plane 
E u mode softens to the point that it also contributes to 
the ground state, producing the r-phase. The frequen- 
cies of the two modes cross at an in-plane lattice con- 
stant roughly in the middle of the r-phase. With further 
increase of the in-plane lattice constant, the A2 U mode 
stiffens to the point that the E u mode is the only mode 
contributing to the ground state, producing the aa-phase. 
The analysis of the Hessian matrix in the ground state 
structures shows that the freezing-in of the polar modes 
eliminates the nonpolar instabilities also seen in Fig. 0] 
most notably the Bi u mode. 

The behavior of pure PT is qualitatively similar. There 
are two unstable polar zone-center modes Am and E Ul 
increasing and decreasing in frequency with increasing 
lattice constant, respectively. The crossing of the two 
modes is likewise very close to the c-aa phase boundary, 
where the r-phase appears when shear strain relaxation 
is allowed. 

For comparison with the unstable modes in the super- 
lattice, in Fig. life) we show the in-plane strain depen- 
dence of the zone-center unstable modes of pure PZ. For 
the unstrained superlattice we have previously shown 6 
that the atomic displacements of the lowest E„ mode, 
mainly confined in the PT layer, are very close to the dis- 
placements of the polar unstable mode of PT bulk. Sim- 
ilarly the atomic displacements of the lowest A2 M mode, 
mainly confined in the PZ layer, are very close to the dis- 
placements of the polar unstable mode of PZ bulk. By 
comparing the phonon frequencies in Fig. 01 we see that 
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FIG. 4: Frequencies of the unstable zone-center phonons of 
the centrosymmetric reference structure as a function of in- 
plane lattice constant for (a) the [001]-(PbTiO 3 )i(PbZrO 3 )i 
superlattice, (b) pure PT, and (c) pure PZ. All modes are 
labeled by their symmetry; solid (open) symbols indicate 
modes with displacements perpendicular (parallel) to the c 
axis. Solid and dashed lines connect the plotted points as a 
guide to the eye, indicating polar and nonpolar character of 
the modes, respectively. Vertical arrows indicate the lattice 
parameter of cubic PT and PZ. 



certain modes in the superlattice can be identified with 
modes in pure strained PT and PZ. The PT-like mode 
of the superlattice, the lowest E u mode, has a strain- 
dependent frequency very similar to that of the lowest E u 
mode of pure PT. The PZ-like mode of the superlattice, 
A2 U mode, compares well with the Pl2 U mode of pure PZ. 
The B2« mode in the superlattice is a PZ-like mode with 
O-displacements mainly confined in the PZ layer and its 
strain-dependent frequency closely corresponds to that of 
the B2 U mode of pure PZ. 
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C. Static dielectric tensor 

The phonon contribution to the static dielectric tensor 
computed as a function of epitaxial strain is presented 
for the superlattice in Fig.JSJa) and (c). In (a), we plot 
the component t ZZl as for thin films this is relevant to 
the dielectric screening of an electric field applied across 
the film. It can be seen that this component diverges on 
both sides of the r-aa phase boundary, and is large in a 
substantial fraction of the r-phase. In (c), we plot the 
inverses of the three eigenvalues of the dielectric tensor; 
thus a value approaching zero corresponds to a diver- 
gence of e. In the c and aa phases, all three eigenvectors 
are uniquely determined by symmetry. In the c-phase 
£3 = £zz is the "collinear" response, i.e. for a field applied 
along the direction of the spontaneous polarization, and 
£1 = £2 = ^xx — £yy In the aa-phase, £3 is the collinear 
component, along [110], £2 is the in-plane component per- 
pendicular to the polarization, and e± = e zz is along the 
direction normal to the layers. Two eigenvalues of the 
dielectric tensor (ei = £2) diverge at the continuous c-r 
transition, and one (£1) at the continuous r-aa transi- 
tion, through there the second noncollinear eigenvalue £2 
is also fairly large. 

In PT (Fig. G2» and (d)), the behavior of the dielec- 
tric constant in the c and aa phases is rather similar to 
that in the superlattice, with a near-divergence in e zz as 
the boundary of the aa-phase is approached. However, 
comparison of Fig. Efa) an d (b) show that the enhance- 
ment of the r-phase in the superlattice gives a large e zz 
over a much larger range of epitaxial strain. 

IV. DISCUSSION 

The most striking features of the phase diagram of 
the epitaxially constrained [001]-(PbTiO 3 )i(PbZrO 3 )i 
superlattice are the width of the r-phase, the separation 
of the minimum energy in-plane lattice constants of the 
c and aa-phases, and the linearity of the energy with 
in-plane strain in the r-phase. These have significant im- 
plications for the coherence of lattice-mismatched thin 
films and for the dielectric response. 

From Fig.|3Ja), it can be seen that the elastic energy of 
the superlattice as a function of in-plane lattice constant 
is not parabolic. Rather, as the in-plane lattice constant 
decreases from 7.73 a.u., the energy rises linearly with a 
small slope to the minimum of the c-phase at 7.50 a.u., 
behaving parabolically only for smaller an. This makes 
the film compliant to a broad range of substrates; the es- 
timated critical thickness of fully coherent films is close to 
constant in this range. For tetragonal PbTiC>3, in com- 
parison, within each phase (c or aa) , the energy function 
is parabolic, so the critical thickness behaves as usual, in- 
creasing inversely to the deviation of the in-plane lattice 
constant from the equilibrium value for that phase. 

Based on the computed energies for the various phases 
of the superlattice in Fig. |3Ja), one can predict the oc- 



currence of structural phase transitions with increasing 
thickness, as the system changes from fully coherent to 
partially relaxed and finally to fully relaxed. For systems 
under compressive stress from the substrate, there should 
be a transition from the c phase to the r phase with in- 
creasing thickness, while for systems under tensile stress, 
a transition from the aa phase to the r phase is expected. 
For pure PT, systems under tensile stress are predicted 
to undergo a transition with increasing thickness from 
the aa-phase through a narrow r-phase into the c-phase. 

The present analysis identifies the most stable super- 
lattice structure with the ten-atom cell shown in Fig. ^ 
The possibility remains that there may be instabilities 
to structural distortions that double or further expand 
the unit cell, in the plane, along the stacking direction, 
or both. This may be relevant for making comparisons 
with experimental observations. 

Certain components of the dielectric tensor diverge as 
the phase boundaries are approached. While the polar- 
ization induced by an electric field along the direction of 
the spontaneous polarization is rather modest, the non- 
collinear responses are of great interest. This polariza- 
tion rotation mechanism, through the coupling of polar- 
ization and strain, has been invoked to explain the large 
piezoelectric response in single- crystal relaxor a 9 ! 10 , and 
to explain the large piezoelectric response computed for 
ordered PZT (with the same structure as the superlattice 
considered here) in Ref.0- The r-phase, with its continu- 
ously variable direction for the spontaneous polarization, 
is expected to have particularly large responses. In par- 
ticular, the zz component of £, the normal polarization 
induced by a field along the normal, will be enhanced by 
a noncollinear contribution throughout the r-phase. 

Finally, we consider the origin of the different behav- 
ior of the two systems. First, we discuss the difference 
in equilibrium lattice constant of the c-phase and the 
aa-phase. In a simple perovskite, such as pure PT, the 
difference arises from coupling between the strain and the 
unstable polar mode. PT has an unusually strong strain 
coupling, leading to a c/a of about 1.06 in the tetrago- 
nal ground state; the c-phase and aa-phase equilibrium 
in-plane lattice constants are 7.32 a.u. and 7.45 a.u., re- 
spectively. Analogous information about the behavior of 
pure PZ can be obtained from the published parameter- 
ization of first-principles results^ with the c-phase and 
aa-phase equilibrium in-plane lattice constants being 7.71 
a.u. and 7.83 a.u., respectively. A simple elastic-energy 
model for the superlattice can be constructed as the sum 
of the elastic energies of the two layers. For each con- 
stituent, the elastic constants of the c and aa phases turn 
out to be similar, with the values for PZ about 70% larger 
than those for PT. In the superlattice, the first-principles 
equilibrium in-plane lattice constant of the c-phase is 0.07 
a.u. less than that of the elastic-energy model, and in the 
aa-phase it is 0.05 a.u. greater, leading to a separation 
of the first-principles lattice constants of the two phases 
about twice as large as in the model. This deviation from 
a simple elastic energy model might be due to the thin- 
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FIG. 5: The phonon contribution to the static dielectric tensor component e zz as a function of the in-plane lattice constant for 
(a) the [001]-(PbTiO3)i(PbZrO3)i superlattice, and (b) pure PT. Inverses of the eigenvalues of the phonon contribution to the 
dielectric tensor as a function of the in-plane lattice constant for (c) the [001]-(PbTiC>3)i (PbZr03)i superlattice, and (d) pure 
PT. The lines are a guide to the eye. 



ness of the individual layers and/or to the importance of 
inter-layer interactions beyond lattice matching. 

The stability of the r-phase is favored by in-plane and 
perpendicular unstable polar modes having both a simi- 
lar strength and positive anharmonic coupling. In a sim- 
ple perovskite, such as pure PT, tetragonal strain splits 
the three-fold degenerate mode to give a dominant unsta- 
ble mode for c/a less than 1, and a dominant in-plane 
E u mode for c/a greater than 1. If the ground state is ob- 
tained by freezing in only the dominant unstable mode, 
we expect a discontinuous change in character at the in- 
plane lattice constant corresponding to the lowest-energy 
cubic structure, as seen in PT. Even if the strengths of 
the two modes are similar, the fourth order anharmonic 
term coupling them must be large enough to favor this 
combination over the higher-symmetry structures (z, as 
in the c-phase, or along [110] as in the aa-phase). Its 
magnitude can vary greatly from material to material, 
being enhanced by coupling to particular types of strain, 
especially shear strain. 

In the superlattice, it seems that separating the two 
components of the polarization in two spatially distinct 
regions helps to stabilize the r-phase structure. Even 
if the strengths of the unstable polar modes in the two 
layer are different, a low-energy structure can be obtained 



by combining the perpendicular mode with polarization 
along z in one layer (e.g. the PZ layer) with a mode 
of different character, with polarization along [110], in 
the other layer (e.g. the PT layer). This idealized pic- 
ture will be somewhat modified by electrostatic consid- 
erations, which will tend to minimize the divergence of 
P by polarizing the PT layer along the z direction. 

The usefulness of this simple picture lies in its applica- 
bility to longer period PT/PZ superlattices; such systems 
are also of interest as they are accessible to experimental 
investigation. 



V. CONCLUSIONS 

The phase diagrams for the [001]-(PbTiO 3 )i(PbZrO 3 )i 
superlattice and pure PT with varying in-plane strain on 
a cubic substrate have been obtained from first-principles 
density-functional total-energy calculations and a stabil- 
ity analysis involving the computation of zone center 
phonons and their coupling to strain. At the smallest 
a||, both systems are tetragonal with pseudo-perovskite 
c/a greater than 1 and P along [001] (the c-phase), and 
at the highest a||, both are pseudo-tetragonal orthorhom- 
bic with pseudo-perovskite c/a less than 1 and P along 
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[110] (the aa-phase). However, at intermediate in-plane 
strain, the transitional behavior is quite different. In the 
superlattice, P rotates away from the normal to give a 
monoclinic structure (the r-phase) , and at a higher value 
of an, completes the rotation into the plane of the lay- 
ers. In the monoclinic phase, the dielectric and piezo- 
electric response is high, the sensitivity to applied fields 
and stresses resulting from the ease with which the polar- 
ization can be changed through rotation. In contrast, in 
pure PT changing epitaxial strain results in a discontin- 
uous transition from P along [001] to P in plane if only 
tetragonal strain is allowed, with full coupling to shear 
strain opening a narrow range of r-phase stability. While 
the characteristic divergence of certain components of the 
dielectric tensor is seen as the phase boundaries are ap- 
proached, there is no significant contribution from po- 
larization rotation. For both systems, we relate the ob- 
served ground states to the instabilities of a relaxed cen- 
trosymmetric structure computed via DFPT. From these 
results we conclude that in the superlattice, the presence 



of two unstable modes confined in different layers, of dif- 
ferent symmetry and comparable strength significantly 
enhances the window for the monoclinic structure, which 
suggests new approaches for designing artificially struc- 
tured materials with high dielectric and piezoelectric re- 
sponse. 
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